arXiv: 1509.0733lv2 [hep-ph] 28 Apr 2016 


Violations of discrete space-time symmetries 
in chiral effective field theory 


J. de Vries^ and Ulf-G. 

^ Institute for Advanced Simulation, Institut fiir Kernphysik, and Jiilich Center for 
Hadron Physics, Forschungszentrum Jiilich, D-52f25 Jiilich, Germany 

^ Helmholtz-Institut fiir Strahlen- und Kernphysik and Bethe Center for Theoretical 
Physics, Universitdt Bonn, D-53115 Bonn, Germany 

^ JARA - Forces and Matter Experiments, Forschungszentrum Jiilich, D-52f25 Jiilich, 

Germany 

JARA - High Performance Computing, Forschungszentrum Jiilich, D-52f25 Jiilich, 

Germany 


Abstract 

We review recent progress in the theoretical description of the violation of discrete space¬ 
time symmetries in hadronic and nuclear systems. We focus on parity-violating and time- 
reversal-conserving interactions which are induced by the Standard Model weak interaction, 
and on parity- and time-reversal-violating interactions which can be caused by a nonzero 
QCD 6 term or by beyond-the-Standard Model physics. We discuss the origins of such inter¬ 
actions and review the development of the chiral effective field theory extension that includes 
discrete symmetry violations. We discuss the construction of symmetry-violating chiral La- 
grangians and nucleon-nucleon potentials and their applications in few-body systems. 



1 Introduction 


Observations of the violation of the discrete space-time symmetries parity (P), charge conju¬ 
gation (C), and time reversal (T), have played an important role in the construction of the 
Standard Model (SM). The observation of P violation (PV) in nuclear /3 decay gave rise to 
the V-A structure of the weak interaction, while observations of CP violation (CPV) lead to 
the prediction of a third generation of quarks. The pattern of discrete symmetry violations is 
nowadays accurately described by the SM. Nevertheless, a number of outstanding issues remain. 

Although flavor-conserving PVTCQ in the SM is well understood at the level of elementary 
fields, its manifestation at the nuclear level is not clear. The SM predicts PVTC interactions 
between quarks, but their appearance at the hadronic level is masked by the nonperturbative 
nature of low-energy QCD Bi- In fact, a consistent framework that satisfactorily describes all 
existing PVTC data does not exist at the moment. The problem of how the PVTC four-quark 
operators appear at the nuclear level thus not only tests the SM weak interaction, but also 
provides an interesting testing ground for nonperturbative QCD [^. 

So far, all observations of CPV in the K and B sector are explained by the phase of the 
quark mixing matrix [^. However, this phase appears to be insufficient to explain the universal 
matter-antimatter asymmetry which requires additional CPV [^. The SM contains a second 
CPV source in the form of the QCD 9 term, which is flavor conserving and severely constrained 
by the non-observation of the neutron electric dipole momenf|^ (EDM) j^. Because the CKM 
mechanism predicts extremely small EDMs [^|^, the search for EDMs is a promising method to 
look for beyond-the-SM (BSM) physics. Current EDM measurements probe high-energy scales 
comparable to the reach of the Large Hadron Collider. To interpret current EDM limits, and 
hopefully future nonzero measurements, a solid understanding is needed of how flavor-diagonal 
P and T violation (PVTV), is manifested at the hadronic and nuclear level [^11|. 

The third possibility of discrete symmetry breaking consists of P-conserving but T-violating 
(PCTV) interactions. Such interactions are only induced in the SM via combinations of PVTC 
and PVTV effects and are therefore very small. Searching for PCTV effects then seems a 
promising way to look for BSM physics. However, any BSM PCTV interactions can mix with 
SM PVTC resulting in PVTV interactions that are strongly constrained by EDM experiments 
|12| - [T4] . Direct PCTV measurements involve scattering processes whose precision is outmatched 
by EDM searches. Nevertheless, scenarios exist where PCTV effects could avoid the EDM 
constraints [^. We do not further discuss PCTV, but refer to Refs. 16-18 for some recent 
studies. 

Traditional approaches of discrete symmetry breaking in nuclear systems are mostly based on 
one-boson-exchange (OBE) models that describe symmetry breaking in terms of various meson- 
nucleon interactions. Although these models provide a useful parametrization, it is difficult to 
link them to the underlying physics and it is unclear whether they capture all aspects of the 
problem. The understanding of low-energy strong interactions has increased tremendously in the 
last decades by the use of effective held theories (EETs), in particular chiral EET (xEFT) (I9p0 


the low-energy EFT of QCD. In this work we review the extension of x^FT that includes the 
breaking of discrete space-time symmetries. 

The xFFT approach has a number of advantages over OBE models. First of all, there exists 
a direct link to the microscopic theory at the quark-gluon level. Another advantage is that the 


^We assume the CPT theorem to hold which implies that CP conservation is equal to T conservation (TC). 
^Nonzero EDMs require the violation of both T (and thus CP) and P. 
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symmetry-conserving and -violating interactions among the relevant degrees of freedom, pions 
and nucleons, are obtained within the same framework which allows for consistent calculations. 
Such calculations can be improved order by order in a controlled expansion within yEFT. Here 
we discuss the PVTC and PVTV potentials up to next-to-next-to-leading order (N^LO) where 
they have a richer form than the corresponding OBE models. Finally, the chiral approach can 
be extended to multi-nucleon or electromagnetic interactions allowing for a unified treatment of 
various different observables. Some of these observables are discussed below. 

In this work we study low-energy PVTC and PVTV interactions in the framework of xEFT. In 
Sec. we discuss the structure of symmetry-breaking interactions at the quark-gluon level both 
from within the SM and beyond. In Sect.we briefly discuss xEFT and describe the construction 
of the PCTC, PVTC, and PVTV nucleon-nucleon (NN) potentials, focusing on their differences 
and similarities. We also discuss the current status of the low-energy constants (LECs) appearing 
in the Lagrangians. In Sect. we focus on several applications of the obtained potentials. In 
the PVTC case, we discuss observables in few-body systems and a recent fit of the weak pion- 
nucleon coupling to existing data. In the PVTV case, we report on calculations of EDMs of light 
nuclei and their possible implications in the search for BSM physics. We summarize and give an 


where more details on the topics addressed below can be found. 


outlook in Sect. Finally, we refer to several recent reviews on PVTC and PVTV 10,21 


2 Symmetry violations at the microscopic level 


The various possibilities of discrete symmetry breaking have very different origins. It is impor¬ 
tant to investigate these origins in order to constrain the operators appearing at lower energies 
in the xEFT Lagrangian. In principle, we could start the analysis directly at the hadronic scale 
by writing down all possible interactions among hadronic degrees of freedom that break the 
discrete symmetries under investigation. The disadvantage of such an approach is that valuable 
information about the hierarchy and structure of the hadronic interactions is lost. The chiral 
symmetry properties of the microscopic PVTC and PVTV operators constrain the form of the 
hadronic operators 22-26 . It is therefore wise to start the analysis at a scale where QCD is 


still perturbative, somewhere above the chiral-symmetry-breaking scale Ay ~ 1 GeV. 


2.1 Origins of PVTC interactions 

Within the SM, PVTC arises from the different gauge-symmetry representations of the chiral 
fermions. As a consequence, only left-handed quarks and leptons are sensitive to the charged cur¬ 
rent weak interaction. Both left- and right-handed fermions interact via the neutral current weak 
interaction, but with different strengths such that P (and C) is still violated. Because the phys¬ 
ical weak gauge bosons, and Z, have masses much larger than the typical hadronic/nuclear 
scale, they can be integrated out and matched to effective PVTC four-quark operators. At a 
scale slightly below Mw, the mass of the W boson, the PVTC operators involving the u and d 
quarks can be conveniently written as 


TpVTC 


Gf 

^/2 




- - sin^ Ow 


— sin^ 6w 


V3a1^3 


3 J 
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Here we have defined = q'j^T^'q, and = qj^q in terms of the quark doublet 

q = (ud)'^. Furthermore, Gp — 1-16 • 10“^ GeV“^ is Fermi’s constant, sin^ ~ 0.23 defines 
the weak mixing angle, and we have set the Cabibbo angle to unity. The dots denote operators 
involving heavier quarks. In particular, operators with strange quarks can have important 
consequences for nuclear PVTC effects. The Lagrangian in Eq. Q needs to be brought to lower 
energies via renormalization group (RG) evolution which dresses the coupling constants with 
0(1) QGD factors and induces operators with different color structure [^,28 . 

The three operators all break P but have different chiral symmetry properties, with the 
first transforming as a scalar, the second as an isovector, and the third as an isotensor [22] . 
Each of the operators thus induces different yEFT Lagrangians. Although all operators in 
Eq. 0 are proportional to Gp ^ there is nevertheless a (small) hierarchy in the sizes of 

the couplings. Due the smallness of sin^ 9w — 0.23, the coupling of the isovector (isotensor) 
operator is suppressed by a factor 5 (2) with respect to the first operator. This suppression 
could potentially lead to smaller LECs arising from the isovector operators, but whether this is 
actually the case depends on, so far, unknown strong matrix elements. In addition, operators 
involving strange quarks with identical chiral symmetry properties, such as do not 

scale with sin^ 6w and could overcome the suppression factors 


22 


29 


2.2 Origins of PVTV interactions 

As discussed briefly in the introduction, the SM has two CPV sources. The phase of the quark 
mass matrix predominantly manifests in flavor-changing interactions and is responsible for the 
observed CPV in K and B decays. Via electroweak corrections flavor-diagonal PVTV interac¬ 
tions are induced, but the suppression factors make them immeasurably small with expected 
experimental accuracies. Eor instance, the neutron EDM arising from the CKM mechanism lies 
roughly six orders of magnitude below the current experimental limit [^|^. 

The second source of CPV in the SM appears in the strong interaction in the form of the 
QGD 6 term 30 . CPV arising from the 0 term is closely connected to the quark masses and it 


is useful to consider them simultaneously. Focusing on the lightest two flavors, the non-kinetic 
part of the QGD Lagrangian can be written as 


^m,9 = - [e^^qiMqp + qpMqp) - 


( 2 ) 


in terms of the gluon field strength G, the diagonal real mass matrix M = rh{l -\- er^) where 
fh = [mu + w-d)/2 and e = {mu — md)/{mu + m^), the strong coupling constant gs, the overal 
phase /?, and the angle 0. The interaction proportional to 0 is a total derivative but contributes 
to the action via instantons [^. For yEFT purposes it is useful to perform an anomalous 
U{1)a transformation to eliminate the 6 term in favor of a complex mass term proportional to 
9 = 6 + 2p. A subsequent non-anomalous SU{2)a rotation can be used to cast the complex mass 
term in isoscalar form 


25,31 




qq + e qr^q — 




9qij^q 


(3) 


where we used 9 1 (the current constraint on the neutron EDM forces 0 < 10 32 ). The 


lack of explanation for this extreme smallness is usually called the strong CP problem. A popular 
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solution is the so-called Peccei-Quinn mechanism which dynamically explains the smallness 
of 6 at the cost of a so far unmeasured new particle, the axion. 

Generic BSM models contain additional CPV phases that can induce larger EDMs than the 
SM [^. Since many BSM variants exist, from the point of view of low-energy experiments it is 
preferable to perform a model-independent analysis. Considering the success of the SM, it is 
likely that any BSM physics appears at a scale considerably higher than the electroweak scale. 
This scale separation makes it possible to treat the SM as the dimension-four and lower part 
of a more general EFT containing higher-dimensional operators. Such operators must conserve 
the all-important SM Lorentz and U{1)y x SU{2)l x SU{3)c gauge symmetries that put strong 
constraints on their form [^. For our purposes, the first relevant operators appear at dimension 
six and are suppressed by two powers of the scale, Acpv, where the additional CPV is assumed 
to appear. BSM CPV can be studied in a model-independent way by adding all CPV dimension- 
six operators , 35 . The great advantage is that it is not necessary to choose a specific SM 


extension, but, if so wanted, the approach can easily be matched to specific models 11 


The advantage of starting the EFT analysis at Aqpv is that we make full use of the SM gauge 
symmetries. The disadvantage is that we must evolve the operators to which is more involved 
than in the PVTC case. The RG evolution involves QCD and electroweak corrections 36 -38 
and heavy SM particles need to be integrated out at their thresholds. A full study has been 
performed for operators involving the lightest two quarks [38], while the evolution of operators 
with heavier quarks is more fragmented, see for instance 39-41 . QCD RG typically mildly 


suppresses the operators when evolved to lower scales and mixes several operators, making it 
more difficult to identify the high-energy origin from EDM experiments. 

Despite these difficulties, the number and form of PVTV interactions around A^^, is very limited 
|24| . Being interested in hadronic and nuclear PVTV we neglect (semi-)leptonic operators. The 
first operators that appear are then the quarlj^ EDMs and chromo-EDMs (CEDMs) 


TpvTV = — 2 


- ^ lido + d3T^)icr^'"l5taqG’l;^^ , 


( 4 ) 


where denotes the photon field strength and ta the generators of SU{3). Although these 
operators appear to be dimension five, SU (2 )p gauge invariance forces them to be coupled to the 
Higgs field at high energies making them effectively dimension six. After electroweak symmetry 
breaking, the Higgs field takes on its vacuum expectation value (vev) resulting in a quark mass 
(the actual value of the mass is model dependent and could be larger or smaller) appearing in 
the scalings do, 3 )C^o ,3 f^^CPV 

The next operator is the Weinberg operator, also called the gluon CEDM (gCEDM) [42||43| 

Tpvtv = , (5) 

which is the only purely gluonic CPV operator at dimension six. Finally, there are several PVTV 
four-quark operators 


Tpvtv = Im S [qq qij^q - + Im S + . .. , (6) 


in terms of 5“ = qr'^q and and the dots denote two similar operators with different 

color structure [^. The coupling constants dw, ImS, and ImH parametrize BSM physics and 

®Here we focus on the lightest two quark flavors and neglect strange q(C)EDMs. 
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scale as A^^py. The operators have very different origins. The operators proportional to dw 
and Im S conserve all gauge symmetries and appear directly at Acpv 44 . Because they are 


invariant under chiral symmetry transformations their low-energy consequences are similar and 
in what follows we only focus on the gCEDM. The operator proportional to ImH is induced 
after electroweak symmetry breaking via a exchang^ between quarks 45 . 


3 Chiral Lagrangians and nucleon-nucleon potentials 

Now that we have identified the structure of PVTC and PVTV interactions slightly above 
A^, we must understand how these interactions manifest themselves at lower energies where 
QCD becomes nonperturbative. To do so, we use the framework of yEFT. By constructing 
the most general Lagrangian in terms of the relevant low-degrees of freedom that incorporates 
the symmetries of the microscopic theory (QCD), we obtain an EFT, called chiral perturbation 
theory (xPT), which is the low-energy equivalent of QCD. The big advantage of xPT is that 
observables can be calculated in perturbation theory. The expansion parameter is p/A^, where 
p is the momentum scale appearing in the process. The nonperturbative nature of low-energy 
QCD is captured by the low-energy constants (LECs) associated with each interaction. These 
LECs need to be htted to data or calculated with nonperturbative methods such as lattice QCD. 
xPT was originally formulated for mesons, but has been extended to the nucleon and multi¬ 
nucleon sectors, where it is usually called yEFT. A big success of xEFT is the derivation of 
the structure and hierarchy of multi-nucleon interactions. The strong NN potential has been 
derived up to N'^LO 46 and successfully describes many few-body processes. Reviews on xPT 


and xEFT can be found in Refs. (19 20 47 


A special role in xEFT is played by the pion triplet. Pions emerge as Goldstone bosons 
of the spontaneously broken chiral symmetry to its isospin subgroup SU{2)l x SU{2)ji —)• 
SU{2)j. Interactions of Goldstone bosons are proportional to their momenta which explains the 
perturbative nature of xPT at low energies. Because chiral symmetry is only an approximate 
symmetry of QCD, being broken by the quark masses and charges, the pions obtain a small 
mass and become pseudo-Goldstone bosons. However, the smallness of the symmetry-breaking 
terms ensures that they can be incorporated in the expansion. The xEFT Lagrangian is then 
obtained by constructing all chiral-invariant interactions and all interactions that break chiral 
symmetry in the same way as the chiral-breaking sources at the quark level. In principle, an 
infinite number of interactions exist, but they can be ordered by the chiral index A = d+nj2 — 2, 
where d counts the number of derivative^ and quark mass insertions (as ~ ~ each 

quark mass insertion increases d by 2) and n the number of nucleon fields [50[|51] . We use 
this definition of the chiral index throughout this review (Weinberg power counting), but refer 
to Ref. 52 for discussions of possible enhancements (i.e. lower chiral indices) for short-range 


nuclear current operators. Because the PVTC and PVTV interactions are associated with very 
small parameters, e.g. G^A^ ~ 10“^ and 6 < they can be included in the expansion as 

well. In fact, the interactions are so weak that only x^FT operators linear in the symmetry¬ 
breaking parameters must be considered. 


"^This exchange involves a left- and right-handed current and the resulting operator is usually called the 
four-quark left-right operator (FQLR). 

^Because tun / i® ^ small number, derivatives acting on nucleon fields are in principle not suppressed. 


In this work we resolve the issue by considering heavy-baryon yPT [48[|49| , where the nucleon is treated nonrela- 
tivistically and the large nucleon mass mjv is removed from the nucleon propagators. 
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An important ingredient in few-body calculations is the NN potential which can be derived 

Weinberg showed |53| that the importance of a connected 


from the xEFT Lagrangians 53,54 


irreducible diagram with N nucleons, L loops, and Vi insertions of an interaction with chiral 
index Aj is given by the index 


u = -A + 2N + 2L + Y,Vi^i ■ ( 7 ) 

i 

The index v is bounded from below which allows for an order-by-order derivation of the PCTC, 
PVTC, and PVTV NN potentials. 


3.1 PCTC chiral interactions 


We begin by listing the most important interactions that conserve P, T, and C and originate in 
the QCD Lagrangian. The construction of the chiral operators is well known and explained in 
various reviews [20||47| and here we only discuss some of the relevant operators. Because chiral 
symmetry is realized in nonlinear fashion in xEFT, each interaction is associated with other 
interactions involving more pions. For simplicity, we only focus here on the operators with the 
least number of pions. The operators with the lowest chiral index A = 0 are given by 

^PCTC ~ ~ ~ 

2 2 Ft, 

-]^Cs{NN){NN) + 2Ct{NS^N){NS^N) +... , (8) 


in terms of the pion triplet if, the nucleon doublet N = (pn)'^, the nucleon velocity and 
spin S^, the pion decay constants Ft, = 92.4 MeV, and three LECs gA, Cs, and Ct- The dots 
denote associated interactions with more pions such as the Weinberg-Tomozawa vertex. Apart 
from the pion mass term, which involves pions without a derivative, all operators originate in 
the chiral-invariant part of the QCD Lagrangian. 

What is important to notice is that chiral symmetry and P and T conservation ensure that the 
leading-order (LO) pion-nucleon vertex comes with a derivative, while the LO NN interactions 
do not. Therefore, both interactions have chiral index A = 0 and contribute to the LO NN 
potential (with u = 0). At NLO, with chiral u = 2 and thus a relative suppression of (p/A^^)^, 
there are additional contact terms with two derivatives, and two-pion-exchange (TPE) diagrams 
involving LO vertices. At N^LO additional TPE diagrams appear which involve yr-Tr-nucleon 
interactions with chiral index A = 1, the so-called c, interactions 


55 , which also contribute to 


three-body forces appearing at the same order. The schematic hierarchy of the PCTC potential 
is shown in the left panel of Fig. We do not discuss higher-order corrections to the Lagrangian 
or potential, but refer to Refs. |46l|56| for the full expressions. 


3.2 PVTC chiral interactions 

The four-quark interactions in Eq. 0 give rise to PVTC hadronic interactions. The three four- 
quark interactions transform, respectively, as a chiral singlet, an isovector (A/ = 1), and an 
isotensor (A/ = 2) [22] . This implies that the first operator can only induce pionic operators 
with derivatives, while the other two can lead to non-derivative pionic interactions. Only the 
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Figure 1: Chiral expansion of the PCTC (left) and PVTC (right) nucleon-nucleon force up to N^LO. Solid 
and dashed lines denote nucleons and pions, respectively. Circles and triangles represent, respectively, 
LO PCTC and PVTC interactions while circled vertices represent NLO interactions. Corrections to the 
one-pion-exchange potential and three-body forces are not shown. 


isovector interaction leads to an 


interaction with chiral indesi^ A = — 1: 


/’(“!) _ 

■^PVTC ~ 


-^A(7f X t)^N 
V2 


(9) 


where h-,^ is the (in)famous weak pion-nucleon coupling constant whose size is still unknown, 
despite significant theoretical and experimental effort. In contrast to the PCTC case, there 
appear no NN contact interactions at LO. Such terms require at least one derivative and have 
chiral index A = 1. Thus, the LO PVTC potential consists only of a one-pion-exchange (OPE) 
contribution: 


w(-l) _ 
•^PVTC ~ 


9 A hn 

2V2F^ 


z(ri X T2f 


(<?i + ^2) ■ q 

n^2 _|_ g2 ’ 


( 10 ) 


in terms of the nucleon spin < 71^2 and the momentum transfer flowing from nucleon (1) to nucleon 
(2): q = p — p' [q = |^|), where p and p' are the relative momenta of the incoming and outgoing 
nucleon pair in the center-of-mass frame. 

Because the LO potential consists of a single term, it might be expected that this term 
dominates hadronic and nuclear PVTC. From the existing data it should then be possible to fix 
the size of h-,^ from which other processes can be predicted. Unfortunately, the situation is more 
complicated for two main reasons: 


• The LO potential changes the total isospin of the interacting nucleon pair and therefore 
does not contribute to PVTC effects in proton-proton (pp) scattering. As a signihcant 
part of the nonzero PVTC measurements has been made in this process, higher-order 
corrections are required to analyze the data. 


• The division of the potential into LO, NLO, .... , is based on an expansion in p/A.^- 
However, as discussed in Sec. |2.1t the isovector four-quark operator is suppressed by a 
factor sin^ 6w ~ 1/5. Large Nc arguments indicate that /itt could even be further sup- 
Thus, formally higher-order corrections might be larger than expected 
because of dimensionless factors not captured by the chiral counting. 

® Although, this PVTC interaction has a lower chiral index than the LO PCTC interactions in Eq. Q, this of 
course does not imply that the interaction is larger. All PVTC interactions are proportional to Gf which greatly 
suppresses their effects. The chiral index should only be used to probe the relative importance with respect to 
interactions with the same transformations under C, P, and T. 


pressed 29,57 
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A great advantage of xEFT is that higher-order corrections can be systematically calculated. 
The NLO potential (with v = 1) was first obtained in Ref. |23| and shown to consist^ of TPE 
diagrams proportional to /itt and NN contact interactions [61| . The TPE contributions suffer 
from ultraviolet divergences that are absorbed, together with the associated scale dependence, 
by the contact terms. We show results using spectral function regularization with cut-off A 5 
(varied between 500 and 700 MeV in the discussions below) which was introduced in Ref. 62 to 


improve the convergence of the PCTC potential. In this way, the PVTC and PCTC potential are 
regularized in the same way. By taking the spectral cut-off Kg —>■ 00 , the results in dimensional 
regularization are retrieved. The TPE contributions are then given by 


23,63,64 




2y/2F^ (dTrP; 


[i(ri X T 2 f{ai -F 0^2) • q] “ l) ^(g, kg) 


2V2F^ (47rF^)2 

in terms of the loop functions 


+ .. 1 [Kf +f 2 f{^i X <72) • q] L{q,kg) , 


( 11 ) 


u = ^/q^+Aml, L{q,kg) = ^ log ’ s = ^k^ - 4ml . (12) 

The first term in Eq. (0 has the same spin-isospin structure as the OPE potential and is 
therefore not very interesting. However, the second term induces ^5o -f-)- ^Pq transitions and 
give the first contributions to pp scattering. 

At the same order as the TPE diagrams, we find fiv^ NN contact interactions. These can 
be written in various ways [61l|65], and here we use the following parametrization (64] 




( 1 ) 


PVTC 


Co 


+ 


F.Al 


((?! - (T 2 ) • (P + P) 


n , ^ (b + ^ 2 )^ , ^ Ti • T2 - 3rfT|\ - n - 

Cl -h C 2 ^-h C 3 -^- 1 i{ai X 0 - 2 ) • O' 


C4 


+ WT2-*('B ^ T 2 ) (cti + a2)-q 


F.Al 


(13) 


All together, the NLO PV potential depends on six LECs which need to be fitted to experiments 
or calculated with nonperturbative techniques. 

At N^LO {u = 2) several additional TPE diagrams appear, which have been calculated in 
Ref. (64| . The first part involves no new LECs and is proportional to tt/itt C 4 and has the same 


spin-isospin properties as the second term in Eq. (11). Because of the large size of C 4 ~ 3.4 
GeV-i ^ 


56 


explained by underlying A and p-meson resonances, and the enhancement by a 
factor TT, this term can be expected to dominate the N^LO potential unless h-,^ is very small. 
The second part depends on five new PVTC pion-nucleon and pion-pion-nucleon LECs [22] , 
which will be difficult to fit to the scarce data. At this order we encounter the first contributions 
to PVTC three-body forces which have not been studiec^ and partially depend on new LECs. 
The hierarchy of the PVTC potential is sketched in the right panel of Eig. 

^NLO corrections to the OPE potential can be absorbed in redefinitions of the LO and NLO LECs |59||60|. 

®This number can be understood by noticing that there are five possible S P couplings. One A ^Pi 
transition, one transition, and three, one for each value of mt, ^So <-4 ®Po transitions. The isospin 

properties of the PVTC four-quark operators are rich enough to induce all five NN operators at the same order. 

®PVTC three-body forces in pionless EFT were studied in Ref. 66 . 
























LEG 

DDH ‘best’ value 

DDH range 

KMW 


0 . 46 - 

10“® 

(0.0 ^ 1.2) • 10“® 

0.10 • 

10“® 

Co 

4 . 7 - 

10“® 

(- 5.0 ^ 13 ) • 10 “® 

0.89 • 

10“® 

Cl 

1.2 • 

10“® 

(- 2.5 ^ 4 . 5 ) • 10 “® 

0.11 • 

10“® 

C2 

-2.2 • 

10“® 

(- 5.0 ^ - 0 . 2 ) • 10 “® 

-0.66 • 

10“® 

C3 

1.0- 

10“® 

(0.8 ^ 1.2) • 10“® 

0.41 • 

10“® 

C4 

0.25 • 

10“® 

(- 0.1 ^ 0 . 7 ) • 10 “® 

- 0.05 • 

10“® 


Table 1 : Predictions 64 of the LECs C, based on resonance saturation and two different sets of DDH 


parameters. The first two columns correspond to, respectively, the best values and range from Ref. 
while the third column is based on the predictions of Refs. 


29 57 


67 


3.2.1 Estimates of the PVTC LECs 


Although xEFT allows the determination of the hierarchy and form of the potential, the LECs 
cannot be obtained from symmetry considerations alone. The PVTC LECs, in particular h-,^, 
have been studied with various techniques with varying precision and sophistication. The sim¬ 
plest estimates are obtained from naive-dimensional analysis (NDA) [^, 68 which predicts 
/i-TT ~ Q ~ O(G'FTn-A^) ~ 10“®. This should be seen as an order-of-magnitude estimate which 
roughly probes the size of PVTC in nuclear systems. 

The most applied approach to hadronic PVTC is the one-meson exchange model (or the 
DDH model, after the authors of Ref. [^), in which PVTC is described by the exchange of a 
single pion, p-, or w-meson. The exchange of a charged pion gives rise to the same potential as 


Eq. (10), while exchanges of the heavier mesons give rise to different structure with shorter range. 
The DDH potential then depends on and six constants associated with heavier mesons. In 
Ref. [67] , /itt was estimated using the quark model and SU{6) symmetry, finding the reasonable 
range 0 < /ijr < 1-2 • 10“® and a “best” value of h-,^ ~ 4.6 • 10“^. In similar fashion, a range and 
best value for the other constants were estimated. 

The DDH parameters were calculated using a soliton description of the nucleon in Refs. 
29,57 , finding ~ 1 • 10“^. Similar small values of h-,^ were found in a first lattice QCPP^ 


calculation 69 and a large-Vc analysis 58 . The smaller values of — 10“' seem to be in 
better agreement with data. The absence of a PVTC signal in 7 -ray emission from leads 
to a strong bound on < 1.3 • 10“^ 70 -72 . On the other hand, the measurement of the Cs 
anapole moment 73 prefers a larger value /ijr — 10“® 74,75 , but suffers from larger nuclear 


uncertainties. This confusing situation can be cleared up by analyzing PVTC signals in few-body 
experiments where the nuclear theory is better under control. 

By use of resonance saturation, the contact LECs Ci can by estimated in terms of the DDH 
parameters [^. After integrating out the p- and ca-meson, the DDH potential (apart from the 
OPE part) collapses into the five contact interactions in Eq. (13) plus terms suppressed by 
powers of (p/m,p^^)^. The estimates of the DDH parameters can then be used to predict the 
sizes of Ci- As the xEFT potential contains explicit TPE contributions that are absent in the 
DDH potential, these must be subtracted to get a sensible comparison . By doing so. Ref. [64] 
obtained the estimates in Table for the LECs Ci based on two sets of DDH parameters. The 
predictions vary by roughly an order of magnitude reflecting the large uncertainty. 


^"Disconnected diagrams were not inclnded nor was the result extrapolated to the physical pion mass. 
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qCEDM(Jo,3) 

qEDM((lo,3) 

gCEDM (dw) 

EQLR (ImS) 

90 

-1 

-1 

3 

1 

1 

9i 

1 

-1 

3 

1 

-1 

92 

3 

1 

3 

2 

1 

A 

0 

0 

4 

2 

-2 

01,2 

1 

1 

5 

1 

3 

do,i 

1 

1 

1 

1 

1 


24 


where d 


Table 2: The chiral indices A = d + n/2 — 2 of various PTVT chiral interactions 
counts the number of derivatives, photon fields, and quark mass insertions and n the number of nucleon 
fields appearing in the operator. The chiral indices should only be used to probe the relative strength 
of the interactions for each PVTV source separately. The LECs arising from different sources cannot be 
compared model-independently due the unknown sizes of the parameters 6, do, 3 ) dw, and ImS. 


3.3 PVTV chiral interactions 


We now turn to the Lagrangians induced by the PVTV sources described in Sect. 2.2 


In the PVTC case, the three four-quark operators in Eq. [^were all proportional to Gp and had 
roughly equal strengths. We therefore could construct one PVTC xFFT Lagrangian with a clear 
hierarchy of interactions. Life is not as simple in the PVTV case. As no nonzero flavor-diagonal 
PVTV measurements have been made, we can only put upper limits on the coupling constants 
6, do, 3 , do, 3 , dwi ImS, and ImH and we know nothing about their relative sizes. By focusing 
on specific BSM scenarios it is possible to calculate the hierarchy of parameters, but in this way 
we are no longer working in a model-independent fashion. Consequently, we must construct 
the PVTV xEFT Lagrangian for each source separately. This has been described in detail in 
Refs. p4 -26 and here we focus on the most important findings. 


Because we construct the xFFT Lagrangian for each PVTV source, the same chiral interaction 
can have different chiral indices depending on the PVTV source. To illustrate this we look at 
the PVTV pion-nucleon interactions. The simultaneous breaking of chiral symmetry, P, and T, 
allows for three different non-derivative vertices: 


£pYTv = <?oVr • ttN + qiNtt^N + g2NT^'K^N , 


(14) 


in terms of three LECs 50 , 1,2 • The 9 term in Eq. (§ breaks chiral symmetry as a complex 
isoscalar quark mass and can only directly induce go which then has chiral index ^Q{go) = — 1 . 
To generate 51 , we require an insertion of the quark mass difference, raising the chiral index to 
^d{gi) = 1- Even that is not enough to induce 52 which requires another insertion of the quark 
mass difference or a photon exchang^^ Ag( 52 ) = 3. On the other hand, the qCEDMs in Eq. Q 
break both chiral and isospin symmetry (if |do| — |c? 3 |) and we obtain: Aj ( 50 ) = Aj ( 51 ) = — 1 , 
while a quark mass insertion is needed for 52 : Aj ( 52 ) = 1- Clearly, different PVTV sources 
induce different xFFT Lagrangians. 

The same game can be played for other PVTV interactions. It turns out that for all sources. 


'^Here we follow Ref. 


56 


and count aem/(47r) as (p/A,,-)'* 
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the operators with the lowest chiral index are given b} 


12 


/IpvTV = QoNt ■ ttN + qiNtt^N + misfATr^TT'^ — 2N{do + d3T^)S^Nv^Ffj^i^ 

+CiNNd^{NS^^N) + C 2 NTN • d^iNrS^^N) , (15) 

in terms of seven LECs [24[|35] . Which of these interactions is actually relevant depends on the 
underlying source. In Table the chiral indices for the different sources are summarized. There 
are a few interesting things to point out: 

• For most sources, §2 appears at higher order than go and/or gi. The one exception is the 
qEDM, for which 50 , 1,2 are suppressed compared to do,i. Thus, LO calculations of PVTV 
observables do not require the inclusion of 52 - 

• Similar to the PVTC case, the fact that PVTV pion-nucleon interactions without deriva¬ 
tives exist while PVTV NN interactions require a derivative, ensures the NN potential is 
dominated by OPE diagrams. The one exception is the gCEDM which conserves chiral 
symmetry. As a consequence, its contributions to 50,1 require a quark mass insertion, such 


that OPE diagrams and PVTV NN interactions appear at the same order 77 


The interactions (io,i describe short-range contributions to the nucleon EDMs. For sources 
such as 9 and qCEDMs, they carry a chiral index which is larger by two than the corre¬ 
sponding index for go. As discussed in the next section, this implies that contributions to 
the nucleon EDMs proportional to 50 (one-loop diagrams) and do,i (tree level) appear at 
the same order 


78 


80 


For all other sources, do,i dominate the nucleon EDMs 


81 


82 


For most sources the three-pion vertex A appears at higher order than go and/or 51 and 
its effects will be minor. For the FQLR, A is relatively enhanced leading to additional 
contributions to the NN potential and a PVTV three-body force |24|. 


With the interactions in Eq. (15) it is straightforward to obtain the PVTV NN potential |83]: 
■9Ado, 


VpVTV = 


2E, 

9A /_ 


in ■ '^2) ^'^^2 ,^^2 ^ (*^1 + ^271 • T2) in - n) ■ q- 




+ igi + A /gi(g)) [(<Ti -k CT2)(ri - T2f + {n - n)in + T2f] 


q 


(p + ml 


(16) 


where (q) describes a momentum-dependent correction to the isospin-breaking potential. It 
arises from a one-loop diagram involving A 


24 


15g\mnmN 


1 + 


and is given by 

1 2q‘^/{Aml) 

3q/ (2m^ 


arctan 


2m, 


(17) 


Explicit calculations in light nuclei 84 show that the dominant part arises from the 5 ^ indepen¬ 
dent part in which case the A dependence can be absorbed in 51 —)• 51 -|- A /gj ( 0 ), which we do 
henceforth. In addition, the A vertex induces a three-body PVTV potential which can influence 
EDMs of nuclei with A > 3. Higher-order contributions to the potential, such as TPE diagrams, 
have been studied in Refs. (M,85 . They are not as interesting as their PVTC counterparts 


because the 50 TPE diagrams have the same spin-isospin structure as the LO OPE, whereas the 
sum of 5 i TPE diagrams vanishes. 


^The simultaneous violation of P, T, and isospin symmetry allows for the presence of pion tadpoles C = /tadvra. 
These terms can be eliminated via field redefinitions at the price of additional PVTV interactions in other 


sectors 24 26 


The Lagrangian in Eq. (151 is interpreted as snbsequent to these redefinitions. 
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3.3.1 Sizes of PVTV LECs and the nucleon EDMs 


In order to calculate PVTV observables such as EDMs in terms of the fundamental couplings at 


the quark level, we must know the sizes of the LECs appearing in Eq. (15). The chiral indices 


in Table provide some information about their relative sizes but this is based on NDA and 
not very precise. Various nonperturbative methods have been applied to calculate the LECs, 
see Ref. for an overview. Nevertheless, for most PVTV sources very little is known which 
hampers the interpretation of EDM experiments. 

We begin with the PVTV pion-nucleon and NN interactions. By far, the most is known for 
the 6 term. As was realized in Ref. [^, chiral symmetry relates the PVTV LECs to PCTC 
LECs originating in the quark masses. For instance, go is proportional to 6mN, the strong part 
of the proton-neutron mass splittingj^ A similar, but less precise, relation exists between gi 
and the strong part of the pion mass splitting p5,84 . With lattice-QCD input for the mass 


splittings 187,88 , the following values are obtained 1841186 : 


go = -(15.5 ± 2.5) • 10“^ 6, gi = (3.4 ± 2) • 10“'* 9 


(18) 


such that \gi/go\ 1/5. Owing to the small value of the proton-neutron mass splitting, this ratio 
of couplings is somewhat larger than expected from Table ^ which predicts \gi/go\ 

|85| . This discrepancy reflects the uncertainty of the NDA estimates in Table With similar 
methods, PVTV couplings between heavier mesons and nucleons can be calculated. For example, 
the PVTV r/-nucleon coupling is proportional to the nucleon sigma term [86] , which is nowadays 
determined with high precision 


89 


Values for Ci ^2 can then be estimated using resonance 
saturation: Ci ~ (—8 • 10“^) 0fm^~and C 2 — ( — 1 • 10“^) 9iw?. 

Much less is known for the dimension-six operators. In case of the qCEDMs, values of 50,1 
have been obtained with QCD sum rules [90], but with 0(100%) uncertainties. For the other 


sources no calculations of 50,1 or Ci ^2 exist and at the moment we cannot do better than NDA 
to estimate their sizes. In some cases, chiral symmetry considerations provide some information 
about the relative sizes of the couplings. For example, for the FQLR the ratio of couplings is 
predicted: 50/51 = <J"iA^/( 8 cim^) ~ 0.01 24. 

Next we study do,i which describe short-range contributions to the nucleon EDMs. These 
terms are renormalized by loop contributions proportional to 50 1 and the total EDMs of the 


neutron, dr, 


and proton, dp, up to NLO are given by 


dn — 


dr — 


do — di — 
do -\- di -\- 


egAdo 

egA 


mZ TTirin \ 




50 


2mN) ’ 

27rm7r 


In^ - 

niN 


51 


■Kin-n: 

2m]sf 


(19) 


where e > 0 is the proton charge, dimensional regularization (MS scheme) is applied, and // de¬ 
notes the renormalization scale. Which of these terms are relevant depends again on the source, 
see Table Using (47riV) ~ A^, we see that for 9 and qCEDMs, both doA and go contribute 
at LO. Inserting Eq. (18) into Eq. (19) and g, = rriN, gives dn — do — di — {2 ■ 10“^®) (9 e cm. 
The current experimental constraint on dn M then limits |0| < 10“^*^, assuming no cancellations 


^^SU{3) symmetry also relates go and mass splittings of octet baryons, but these relations suffer from large 
SPlSl-breaking corrections 


86 


The lengthier SU{3) yPl' expressions can be found in Refs. [78 80 92]. 
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between do — di and the loop piece. A more precise constraint requires nonperturbative infor¬ 
mation about the short-range terms. For all other sources, the nucleon EDMs are dominated by 
(io,i and chiral symmetry provides little information about their sizes. 

Lattice QCD can directly determine the LECs (io,i ia terms of the PVTV quark-gluon oper¬ 
ators. In recent years, a lot of effort has gone into determining the nucleon EDMs from the 6 
term 93 ,^. The simulations take place at nonphysical quark masses and in a finite volume 


and xPT expressions are needed to extrapolate to the physical point and infinite volume 95,96 


Based on unpublished lattice data from Shintani et al, Ref. 961971 extracted the following values 


dn = -(2.7 ± 1.2) 


10 ^®ecm 


dp = (2.1 ± 1.2) • IQ-^'^ecm 


( 20 ) 


Ref. 


98 


performed a lattice-QCD calculation of dn by analytically continuing 6 into the complex 


plane. In this way, the Euclidean action becomes real and standard stochastic methods can be 
applied. The following value was found 


dn = -(3.8 ± 1.0) • 10"^®ecm 


( 21 ) 


Finally, Ref. performed a quenched calculation of dn and dp based on the gradient flow for 


gauge fields 100 , finding results consistent with Eqs. (|20|) and (|21|). go was extracted from the 

finding a value a 


calculated electric dipole radius by a comparison to xPT predictions |80, 91 
few times larger than Eq. (18). Considering the limitations of the calculation (quenched and 


800 MeV pion mass), not much can be said about this discrepancy. 

The other well-studied operator is the qEDM. As can be glimpsed from Table the only 
relevant LECs are doq as all other interactions are suppressed by cJem/dTr. The nucleon EDMs 
can be related to the nucleon tensor charges and a recent lattice QCD calculation [101] found 

dn = -(0.23 ± 0m)du + (0.77 ± 0m)dd , dp = (0.77 ± 0.07)d„ - (0.23 ± 0m)dd . (22) 

Much less is known for the remaining dimension-six sources. For the qCEDMs, the nu¬ 
cleon EDMs have been calculated with QCD sum rules with 50% accuracy 102,103 , while 


for the gCEDM and four-quark operators only estimates at the order-of-magnitude level ex¬ 
ist [^ 81 82,104 . Improvements on the hadronic matrix elements could have significant impact 


on falsifying or identifying specific BSM scenarios from EDM measurements, see Refs. Eumi 
for discussions. 


4 Discrete symmetry breaking in few-body systems 

Having discussed the discrete-symmetry-breaking NN potentials, we now turn to observables in 
few-body systems. We mainly focus on systems with two nucleons as these capture the essential 
ideas and avoid the technical difficulties associated with larger systems. The starting point of 
the calculations is the PCTC NN strong interaction from which nuclear wave functions are 
calculated. In most of the literature on discrete symmetry breaking, phenomenological high- 
quality potentials such as Argonne Auig [l06l or Nijmegenll [107] are used which accurately 
describe the NN experimental data-base. The obtained wave-functions are then combined with 
PVTC/PVTV OBE models (such as the DDH model) or with the x^FT potentials. The latter 
method is usually called the hybrid approach. 

Calculations where also the wave functions are obtained from xEFT have been performed 
in various frameworks. At very low energies, the pion can be integrated out. In the resulting 
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pionless EFT |19| , all interactions are described by NN contact vertices which allows for much 
simpler calculations. The absence of pions implies that the EFT only converges in the low-energy 
region E ~ M^/(2m7v) — lOMeV. By integrating out the pion, the hierarchy of forces due to 
the chiral-symmetry properties of the PVTC and PVTV quark-gluon operators is obscured. The 


pionless approach has been applied to several PVTC observables 65 66 108 -110 and is reviewed 
in Ref. [^. 

Another EFT approach is based on the KSW (after the authors of Ref. (nil) power count¬ 
ing. Within this framework, the LO PCTC NN contact interactions are resummed to obtain 
scattering lengths and bound-state properties. Pions are treated in perturbation theory which 
has the advantage that analytic calculations become possible and renormalization is explicit. 
Several PVTC and PVTV observables in the two-body sector have been calculated in this 
way 77,112,113 . Unfortunately, higher-order calculations m show that the KSW expansion 
already fails to converge at fairly low momenta q ~ such that observables associated with 
larger momentum transfer or denser nuclei cannot be treated in this fashion. 

Here we mostly discuss the approach, which we will call yEFT, where pions are treated ex¬ 
plicitly and nonperturbatively and both PCTC and symmetry-breaking potentials are obtained 
from xEFT. Typically numerical differences between hybrid and xEFT calculations are small 
if observables are mainly determined by long-range physics (i.e. pions). Although the PCTC 
xEFT potential can be calculated in perturbation theory, the potential itself must be iterated 
to all orders to calculate NN scattering lengths, phase shifts, and bound state properties. Tech¬ 
nically, the task is to solve the non-relativistic Lippmann-Schwinger (LS) equation 


T = V + VGqT , 


(23) 


written here in short-hand notation. T denotes the scattering matrix, Gq the non-relativistic 
propagator of the free nucleon pair, and V the PCTC strong potential. Once T is obtained 
numerically, the effects of the PVTC and PVTV potentials can be calculated in perturbation 


theory. In practice, it can be easier to solve Eq. (23) directly with V the sum of the strong and 
symmetry-breaking potentials. Because the latter are proportional to very small parameters, 
their iteration leads to numerically identical results as first-order perturbation theory. Details 


on the solution of the LS equation in the presence of P violation can be found in Refs. [TT^[IT6 


The main consequence of the PVTC and PVTV potentials is that transitions between states 
with even and odd orbital angular momentum become possible. For instance, the PVTC OPE 
potential in Eq. (10) leads to transitions normally forbidden by P conservation. 


In general, the momentum integral in the LS equation is divergent and needs to be regulated. 
Different choices are of course possible and most results below are based on Ref. (5^ which 
regularized the LS equation by multiplying the total potential by the functior 


15 


V{p,p) exp[—p' ”'/A^"’] V{p',p) exp[—p ”'/A^”] 


(24) 


where n = 3 and A is a momentum cut-off. Recently, a regularization scheme in coordinate space 
has been developed which better preserves the long-range part of the pion-exchange terms 118 


This scheme has not been applied yet to PVTC or PVTV calculations. Much has been written 
about the size of the momentum cut-off, A, (IT^ |120| , which we cannot discuss in detail here. 
Although A can in principle be any high-energy scale, for practical purposes it is optimal to pick 


Ref. |ll7| applies a similar function but uses n = 2 for NLO and higher-order corrections. The PVTC ^EFT 


analysis of Ref. 60 uses a different regulator for the PVTC potential. 
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A of similar size as [56| 117 . The x^FT calculations discussed below all varied A between 
450 and 600 MeV to get an estimate of (part of the) theoretical uncertainties. 


4.1 Few-body PVTC processes 

We now turn to the calculations of PVTC observables. For an overview of PVTC experiments 
we refer to Ref. . We focus on the observables in the two-body systenj^ for which nonzero 
signals have been measured. Afterwards, we briefly discuss more complicated systems. 


4.1.1 Proton-proton scattering 

The first observable we discuss is the pp longitudinal analyzing power (LAP) which vanishes if P 
is conserved. It is defined as the difference in cross section of scattering between an unpolarized 
target and a beam of positive and negative helicity, normalized to the sum of cross sections. 
The existing experiments measnred the LAP over a certain angular range (from 9i to 62 ) and 
report the integrated asymmetry: 


Al(F,0i,02) 


dQ.{da^ — da-) 
dQ{da+ + da-) 


(25) 


Al has been measured for beam energies of 13.6 MeV [123] , 45 MeV [124 , and 221 MeV 125 


The experiments at 13.6 and 45 MeV are scattering experiments which measured the pp LAP 
over an angular range of, respectively, 20°-78° and 23°-52° (lab coordinates) 


Al{13.6 MeV) = (-0.93 ± 0.21) • 10“'^ , Al( 45 MeV) = (-1.50 ± 0.22) • 10“^ . (26) 


The experiment at 221 MeV is a transmission experiment and effectively measures the LAP over 
the whole angular range modulo a small opening angle [ID 

AL(221MeV) = (0.84 ± 0.34) • 10“'^ . (27) 


Traditionally, it was taken that the pp LAP does not depend on the weak pion-nucleon 
coupling /ijr 115,126 as the OPE potential does not contribute to pp scattering. In yEFT, the 


first non-vanishing contributions appear at NLO and consist of TPE contributions and short- 
range NN interactions. The former depend on /i^r and the latter on the combination of LECs 
C = —Co + Cl + C 2 — Co that gives rise to ^Sq -h- transitions. 

The calcnlation of Al in terms of and C is complicated by the Coulomb interaction. This 
problem has been carefully studied in Refs. [115 ,116 . The low-energy scattering experiments do 
not measure over small angles and the Coulomb interactions plays a minor role. The transmission 
experiment is more affected. The energy of 221 MeV was chosen because the contribution from 
j = 0 transitions are proportional to the sum of the strong phase shifts Ji which vanishes 

around 220 MeV 126ll27 . Therefore ^2,(221 MeV) is sensitive to higher partial waves (mostly 


^P 2 ^D 2 transitions that are proportional to h-jr) and is complementary to the low-energy 

experiments. However, the Coulomb interactions shift the vanishing of the j = 0 phase shifts to 
a somewhat larger energy [TT^ . Combined with the larger uncertainty of the yEFT potentials 
at higher energies, the data at 221 MeV unfortunately does not pin down a precise value for h.,^. 


“We do not discuss observables in the one-body sector such as the proton anapole moment 


for on-shell photons but provides a radiative correction to PVTC electron-proton scattering 122 . 


121 


that vanished 
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The LECs were fittecf^ in Refs, 
from Refs 


56 


116 


and 60 using, respectively, N^LO strong potentials 


and Ea- Ref. |116[ found the strongly correlated (see Fig. values 
K = (1-1 ± 2) • 10"® , C = (-6.5 ± 8) • 10"® , 


116 . Ref. 64 


(28) 

in excellent agreement with Ref. [60] . The uncertainty is almost completely determined by the 
lack of data and the experimental uncertainties, while uncertainties due to cut-off variations are 
much smaller. An additional experiment around 125 MeV would be very beneficial in reducing 

investigated the contributions to from the N^LO 
|3.2[ unless is very small, the dominant part of the 
N^LO potential is expected to be proportional to Including this correction, the fit slightly 

improves h-,^ = (0.8 ± 1.5) • 10“® and C = (—5.5 ± 7) • 10“®, but does not strongly affect the 
values of the LECs indicating that the yEFT expansion is converging satisfactorily. The other 
parts of the N^LO potential depend on additional PVTC LECs and more data is needed before 
they can be htted. 


the uncertainties of the fit 
PVTC potential. As discussed in Sec 


4.1.2 Neutron capture on the proton 


We now discuss the LAP, a.y, for the process^® np —)■ dy which is defined as 


da+{6) - da-{9) 

^ da + (9) + da-{9) ^ 


(29) 


with da±{9) the differential cross section for neutrons with positive/negative helicity and 9 the 
angle between the photon momentum and the neutron spin. For a long time only upper bounds 
existed for a.y [129 130 , but recently a preliminary result was reported [131 


= (-7.14 ±4.4) • lO" 


(30) 


only two standard deviations away from zero. The result is based on a subset of the data taken 
by the NPDGamma collaboration. A full result with uncertainty ~ 10"® is expected soon. 

The great interest in measuring is that, in contrast to A^, it does depend on the LO 
PVTC potential and is therefore more sensitive to h-jt- It has been studied in detail in the DDH 
framework [132 , hybrid calculations [133 ,134 , pionless EFT 108 , KSW-counting [111 , and re¬ 
cently yEFT 135 . A new aspect of is that it depends on PCTC and PVTC electromagnetic 


currents, which can be calculated along the same lines as the NN potentials. The experiment 
takes place at thermal energies where the capture cross section is dominated by the nucleon 
isovector magnetic moment, while NLO two-body currents contribute at the 5% level. In com¬ 
bination, with the N®LO yEFT potential 56 , these currents predict a capture cross section 

The remaining 4% dis- 
These have not been 


319 ± 5 mb 135 , compared to the experimental 334.2 ± 0.5 mb 136 


crepancy should be due to higher-order effects such as TPE currents 137 


included in the analysis as the corresponding PVTC TPE currents have not been calculated. 

Apart from the above ingredients, the LO calculation of requires the nucleon convection 
current (from gauging the nucleon kinetic energy) and PVTC OPE currents proportional to 


115 


The analogous analysis of DDH parameters was performed in Ref. 

The inverse process yd —> np could provide complementary information, see for instance Refs. 


110 128 


^®Note the inconsistenty in the chiral order of the potential and currents as is typical in these kind of calculations. 
Consequently, the uncertainty is dominated by missing higher-order currents. 
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Figure 2: The allowed ranges for the LECs and C (in units of 10“®). The black ellipses correspond 
to fits to the combined pp and np data with total = {2, 3,4}. The red (dashed) ellipses are the same, 
but using the expected future experimental uncertainty (at 10“® level) for the ftp —>■ dj measurement, 
while keeping the central value of Eq. (30). The shaded regions shows the constraints on from 
data and a lattice-QCD calculation of 


^7 then arises from an interference between the isovector magnetic moment and: 1) the 
convection current and the LO PVTC potential, 2) the PCTC OPE currents and the LO PVTC 
potential, 3) the PVTC OPE currents. The individual contributions only have a minor depen¬ 
dence on the cut-off used to regulate the LS equation, but the sum suffers from a relatively 
larger dependence due to cancellations [133, 135 . Ref. [135 


found 


a.^ =-(0.11 ±0.05) . (31) 

The central value agrees with DDH and hybrid results [132 134 that applied the Siegert theo¬ 
rem to relate the electric dipole currents to the one-body charge destiny. The observed cutoff 
dependence is likely to be reduced by the inclusion of higher-order corrections to the PCTC 
and PVTC exchange currents Although these can be partially taken into account by the 


Siegert theorem, this approach most likely overestimates the accuracy of the calculation as not 
all exchange currents are taken into account. 

The analysis oi Ap and can now be combined to extract more precise values of h.,^ and C. 
Contours of total = 2, 3,4 (xmin — 0.7) in the — C plane are plottecj^ in Pig. Varying 


the cut-offs and considering the whole allowed range. Ref. 135] conservatively extracted 
K = (1.1 ± TO) • 10"® , C = (-6.5 ± 4.5) • 10“® 


(32) 


The fits indicate that small values of /itt ~ 10 ^ are still consistent with the data, but that larges 
values of /ijr ~ (5 — 10) • 10“'^ are preferred. Such values are in conflict with the bound on 

/i,r turns out to be very small, the OPE and TPE currents are suppressed and the dominant contribution 
a~, would depend only the PVTC contact LEG C 4 : a~f = —(0.055 ± 0.025) Ca ]135| . A contribution at the same 
order from a PVTC pion-nucleon-photon interaction is neRlig ible 
^'^The plot is based on based on the N^LO yEFT potential 
and As = 600 MeV. 


56 


134 


with intermediate cut-off values: A = 550 MeV 
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gamma-ray emission, h-,^ < 1.3 • 10“'^, and lattice and model calculations of ~ 10“'^. The 
upcoming increase in sensitivity of the measurement will tell whether small values of are 
consistent with few-body experiments as can be seen from the red (dotted) ellipses. 


4.1.3 Additional PVTC observables 


A fairly large set of few-body PVTC observables could provide additional input on nuclear 
PVTC [^. Examples in two- and three-body cases are the spin rotation in np or nd and the 
LAPs of fid scattering and nd which have been discussed in DDH |l28yi38], pionless 


EFT [l09| , hybrid calculations 134 
the spin rotation angles gives 


and, recently in yEFT 60 . The yEFT analysij^ of 


-^{np) = (1.31 ± 0.05)/i^ (0.20 ± 0.01)00 - (0.23 ± 0.01)01 - (0.44 ± 0.01)03 - (0.09 ± 0.01)04 


^(nd) = (2.20 ± 0.02)/i^ - (0.08 ±0.01)Oo - (0.19 ± 0 . 01)04 , 

in units of rad/m. Another interesting observable is the LAP of the reaction n ± ^He —?■ p ± 

Al = -(0.14±0.01)/i^±(0.017±0.003)0o-(0.007±0.001)0i±(0.008±0.001)02±(0.018±0.002)04 

(34) 

These result s are based on the PCTC potential of Ref. [ll7 supplemented by N^LO three- 
body forces [139| . The dependence is dominated by the OPE potential, with TPE contribu¬ 
tions entering at the 10% — 30% level as expected from power counting. The spin rotation angles 
have not been measured yet, but experiments are being considered. The LAP Al is planned to 
be measured with sensitivity 


1.6 • 10"® 140 


Finally, we mention the LAP in pa scattering for which a nonzero value has been measured 
Ai(46MeV) = —(3.3 ± 0.9) • 10“^ 141 . No yEFT calculation exists for this process, which 


is unfortunate as the measurement could have important implications for the extraction of the 
PVTC LECs. Ref. [142 found Al = —0.34/i,r ± • ■ •, based on an optical model and the dots 


denote contributions from other DDH parameters. The difficulty in the yEFT calculation lies 
in the five-body problem and the presence of Coulomb interactions. The great progress in ab 
initio few-body calculations in recent years, see for instance the calculation of aa scattering with 
nuclear lattice EFT 143 , might make it possible to analyze the ffa LAP with similar methods. 


The same applies for the na spin rotation [144 . 


(33) 


4.2 PVTV effects in few-nucleon systems. 

The main motivation for the study of PVTV in few-nucleon systems are the plans to measure 
EDMs of light nuclei in storage rings. Traditional EDM experiments essentially consist in looking 
for a change in the spin precession of the system in the presence of electromagnetic fields. A 
charged particle at rest in an electric field would quickly escape the experimental apparatus. 
This is not true for charged particles confined in a storage ring and it was realized that EDMs 
of charged particles could be measured in such a setup [145 . The JEDI collaboration |146 147 


^Ref. |60| uses a different parametrization of the PVTC NN LECs. 


been translated to Eq. (|13|) via the relations: CY"' = Co/2 — C i, C Y''' = Co/2 


CY''' = C 3 I 2 , where C; denotes the LECs appearing in Ref. 


60 


The results in Eqs. (331 and (341 have 
(yViv ^ ^ _C 2 / 2 , and 
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plans to measure EDMs of light ionq^ such as the proton, deuteron, and helion (^He nucleus) 
in a dedicated storage ring with an accuracy exceeding current neutron EDM experiments. 

The advantage of going to systems with more than one nucleon, is that their EDMs become 


sensitive to the PVTV potential in Eq. (16). The nucleon EDMs depend on the PVTV pion- 
nucleon LECs only via loop diagrams associated with the typical m^/( 47 rEV)^ suppression factor. 
The EDMs of nuclei, however, depend already at tree level on the PVTV pion-nucleon LECs 
such that the EDMs of the bound states can be considerably larger than the constituent nucleon 
EDMs [W 148, 149| . As discussed below, this argument only holds for some of the PVTV sources 
at the quark-gluon level such that combined measurements of EDMs of nucleons and light nuclei 
would point towards the underlying PVTV source 35,77. The second advantage of light ions 


is that their EDMs can be calculated to high precision. Calculations of heavier systems such as 
i99Hg and ^^^Ra rely on nuclear models and suffer from larger uncertainties. These interesting 
observables are outside the reach of current yEET methods and we refer to Ref. for a 
detailed discussion. 

PVTV effects in processes such as np [150 and nd EH] scattering can be studied along 
analogous lines. However, these effects are most likely too small to be measurable. For instance. 
Ref. [150 concluded that with state-of-the art experimental accuracies, measuring the PVTV np 
spin rotation would fall short of the current neutron EDM sensitivities by roughly three orders 
of magnitude. 


4.2.1 The EDM of the deuteron and helion 

The prospect of measuring EDMs of light nuclei in storage rings, has led to a number of investi¬ 
gations of the deuteron EDM 35,77,85 152, 153j . From the theory point of view, the deuteron is 
interesting because its spin-isospin properties ensure that its EDM has rather distinctive prop¬ 


erties. Most calculations are based on the three PVTV pion-nucleon LECs in Eq. (14) and 


phenomenological models of the PCTC NN potential. A calculations in yEFT using N^LO 
PCTC potentials obtained 

do = (0.94 ± 0m){dn + dp) + (0.18 ± 0.02) 51 efm , (35) 


in good agreement with model and hybrid calculations 35, 153j and a calculation based on 


KSW counting 77 . The independence on go and Ci ^2 can be understood from the fact that 


the deuteron is mostly in a ^Si state. After an insertion of the terms in Eq. (16) proportional 


to go or Cl^ 2 ) the wave function obtains a small ^Pi component. The dominant current arises 
from a coupling to the proton charge and is spin-independent. It therefore cannot return the 
wave function to its ground state and the contributions vanish. A recent calculation of 
the ®Li EDM based on a cluster model found a similar effect and a twice larger dependence 
on gi 154 . Refs. (35,85 investigated higher-order contributions to do proportional to go from 


isospin breaking and higher-order currents, but these were found to be negligible. 

The EDMs of ^He and have been calculated in Refs. [^[^|l55l|l56|. Ref. 84 used the 


^^Schiff’s theorem implies that the nuclear EDM in an atomic system is screened by the surrounding electrons. 
The screening is very effective for light atoms, hence the need to measure the EDMs of light nuclei directly. 
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same xEFT potentials as used for do (supplemented by chiral three-body forces) and founc^*^ 


dsHe = (0-90 ± 0.01) dn - (0.03 ± 0.01) dp 

[-(0.017 ± 0.006) A (0.11 ± 0.01) go -h (0.14 ± 0.02) gi] e fm 
[-(0.04 ± 0.02)Ci-h (0.09 ±0.02)C2] efm-^ . 


(36) 


Apart from a smaller contribution from dp and a tiny dependence on A via the PVTV three-body 
forc^^ the most important difference with respect to do is that there is no isospin selection such 
that go and C \^2 contribute. It is this observation that makes a dsHe measurement complementary 
to that ot do- The dependence on Ci 2 in Eq. (36) is larger than found in Ref. 35 , possibly due to 


the pronounced short-range repulsion of the Auig potential used in that work. The uncertainties 


given in Eqs. (35) and (36) are associated to the nuclear theory and are significantly smaller 


than uncertainties associated to the hadronic theory, i.e., the sizes of the LECs in terms of the 
PVTV quark-gluon operators. This is a big advantage over heavier systems where the nuclear 
uncertainty can be the limiting factor. 

The size of d^ and dsHe with respect to dn and dp depends on the PVTV source. Because the 
nucleon EDMs induced by the 9 term are mostly isovector, see Eq. (20), the sum of 0.94(d„ -|- 
dp) = —(0.6 ± 1.6) • lO“^®0ecm appearing in d^ is uncertain. The two-body contribution is of 
comparable size and could potentially cancel against the one-body terms. We therefore focus on 


the two-body contributions and by inserting the values of the LECs in Sect. 3.3.1, we obtain |84 


do - {0M ±0.01) {dn +dp) = 
dsHe - (0-90 ± 0.01) dn ± (0.03 ± 0.01) dp = 


(0.9 ±0.3) • 10“^® 0 6 cm 
-(1.0 ±0.4) • lO"^®0ecm 


(37) 

(38) 


which provides a clean way to extract 9 from measurements of dn, dp, do, or dsHe- The short- 
range operators Ci ^2 contribute at the 10% level as expected from power counting. With lattice 
calculations of dn and dp, these relation can be used to test for the existence of the 9 term given 
measurements of dn or dp in combination with do 01 dsHe- 

Eor the qCEDM, do and dsHe are expected to be dominated by the 50,1 terms. Due to the 
lack of knowledge about the sizes of the LECs, the exact enhancement oido/{dn ± dp) is unclear. 
QCD sum rules and NDA predict dn ± dp (dn) to be roughly 10% (30%) of the pion-exchange 
contribution to do (dsHe); but more precise statements require lattice input of the PVTV LECs. 

The story is similar for the EQLR but now go is suppressed. In this case, the constituent 
nucleon EDMs are expected to enter at the 10% level such that the ratio dsHe/^D — 0.8 can be 
predicted. Such a signal would be a tell-tale sign of existence of the EQLR which is induced in 
left-right symmetric extensions of the SM (see Ref. [I57l and references therein). 

In case of the qEDM, the situation is simple and the light-nuclear EDMs are dominated by 
the nucleon EDMs. That is, do — dn± dp and dsHe — dn- 

Einally, in case of the gCEDM the situation is most complicated. Power counting indicates 
that all contributions to do and dsfje appear at the same order (apart from A). Explicit 
calculations find a somewhat smaller dependence on ^ 0,1 and Ci ^2 then expected 35,84 , but 
nevertheless it is hard to make predictions. Eor instance, even the relative sign of dn and dp is 
unknown which strongly impacts the interpretation ot djo- Lattice input is direly needed. 


do not give the trit on EDM results as it most likely will not be measured due to its radioactive nature. 
^®Power counting 24|[84 predicted a larger dependence on A and it is unclear why the explicit calculation gives 
a small value. 
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The above discussion shows that measurements of EDMs of light nuclei could isolate the 
underlying PVTV source. A different question is whether this information can be used to learn 
something about the possible SM extension at high energies. In Ref. [11| , various popular 
scenarios of BSM physics, such as two-Higgs-doublet, left-right symmetric, and supersymmetric 
models were investigated in the context of EDMs. It was argued that measurements of the 
EDMs of a few systems could not only distinguish such models from a SM 9 term, but also 
partially separate the models based on the EDM hierarchy they induce. The analysis could be 
significantly improved with lattice calculations of the PVTV LECs. 

5 Final remarks 

We have reviewed the breaking of discrete space-time symmetries in strongly-inter acting systems. 
We have focused on flavor-diagonal PVTC and PVTV interactions which have very different 
origin and experimental signatures. The SM weak interaction induces flavor-diagonal PVTC 
four-quark interactions whose structure are well understood. Elavor-diagonal PVTV interactions 
appear only in the SM in the form of the QCD 9 term which is strongly constrained by neutron 
EDM measurements. This smallness leaves room for PVTV effects from BSM physics which at 
low energies can be parametrized by various interactions of dimension six. Although in both 
cases the form of the discrete-symmetry-breaking operators at the quark-gluon level is clear, 
their manifestation at low energies is obscured by nonperturbative QCD. 

To overcome the problem of low-energy QCD, we have focused on the application of y;EFT 
which allows for the systematic derivation of interaction among the relevant low-energy degrees 
of freedom: pions, nucleons, and photons. The universality of this approach is reflected by the 
fact that both PVTC and PVTV chiral Lagrangians and NN potentials can be constructed by 
essentially equivalent methods. The symmetry-breaking potentials can nowadays be combined 
with PCTC xFET potentials to consistently calculate PVTC and PVTV observables in few- 
nucleon systems. 

In the PVTC case, the LO potential consists of a single interaction proportional to the weak 
pion-nucleon coupling, whose size has been an outstanding issue for a long time. The existing 
data on PVTC in pp scattering and ftp —)• dy allow for an extraction h-,^ = (1.1 ± 1.0) • 10“®. This 
is unfortunately not precise enough to rule out or identify the small values of h-,^ < 10“'^ that are 
calculated in lattice QCD and preferred by experiments on ^®F. The upcoming data on ftp —)■ dy 
will hopefully resolve this issue. In any case, additional measurements combined with consistent 
calculations are needed to confirm the size of h-,^ and fix the values of the LECs appearing in 
the PVTC potential. The goal is to finally get a consistent picture of PVTC nuclear forces. 

In the PVTV case, the hierarchy of the NN potential crucially depends on the underlying 
PVTV source. For sources that break chiral symmetry, such as the 9 term or BSM sources like 
quark chromo-EDMs and the FQLR, the potential is dominated by OPE similar to the PVTC 
case. The sources can be differentiated by the relative sizes of the PVTV pion-nucleon LECs. 
The three chiral-breaking sources predict, respectively, \go/gi\ — 5, \go/gi\ — 1, \go/gi\ — 0.01. 
We have discussed how measurements of the EDMs of light nuclei can be used to probe these 
ratios and disentangle the sources. For other sources, such as the gCEDM, the PVTV potential 
also depends at LO on short-range NN interactions which makes the analysis more complicated. 
Considering the expected reach of future EDM experiments these results can play an important 
role in constraining or finding BSM physics. 
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An outstanding problem are calculations of the symmetry-breaking LECs using nonpertur- 
bative methods such as lattice QCD. Both in the PVTC and PVTV sector, very little is known 
about the sizes of the LECs which hampers the interpretation of the experimental data. We 
have discussed recent progress in calculations of the nucleon EDM arising from the 6 term and 
quark EDMs, but stress that similar calculations for the other PVTV sources would be very 
important. The same can be said for calculations of the PVTC and PVTV pion-nucleon LECs. 

In this work we have focused on few-nucleon systems where the scattering equations can be 
solved with high precision. Many experiments have been performed on heavier systems in which 
PVTC and PVTV effects can be significantly enhanced, see discussions in Refs. i [MUM 


The great challenge for the future is to extend the xEFT framework discussed here beyond the 
few-body regime. Great progress has been made in the last couple of years in performing ab 
initio calculations of medium-heavy nuclei based on chiral PCTC NN interactions 159 160 


It would be extremely interesting to extend these calculations to include symmetry violations. 
As discussed above, an interesting intermediate step would be the ab initio calculation of the 
PVTC pa analyzing power for which a nonzero measurement has been reported. 


Acknowledgements 

We thank all our many collaborators for sharing their insights into the topics discussed here. 
This work is supported in part by the DEG and the NSEG through funds provided to the 
Sino-German GRG 110 “Symmetries and the Emergence of Structure in QGD” (Grant No. 
11261130311). The work of UGM was supported in part by the Ghinese Academy of Sciences 
CAS President’s International Fellowship Initiative (PIFI) grant no. 2015VMA076. 


References 

[1] M. J. Ramsey-Musolf and S. A. Page, Ann. Rev. Nucl. Part. Sci. 56, 1 (2006), hep- 
ph/0601127. 

[2] W. C. Haxton and B. R. Holstein, Prog. Part. Nucl. Phys. 71, 185 (2013), 1303.4132. 

[3] M. R. Schindler and R. P. Springer, Prog. Part. Nucl. Phys. 72, 1 (2013), 1305.4190. 

[4] M. Kobayashi and T. Maskawa, Prog. Theor. Phys. 49, 652 (1973). 

[5] L. Canetti, M. Drewes, and M. Shaposhnikov, New J. Phys. 14, 095012 (2012), 1204.4186. 

[6] C. Baker et al., Phys.Rev.Lett. 97, 131801 (2006), hep-ex/0602020. 

[7] 1. B. Khriplovich and A. R. Zhitnitsky, Phys. Lett. B109, 490 (1982). 

[8] C.-Y. Seng, Phys.Rev. C91, 025502 (2015), 1411.1476. 

[9] M. Pospelov and A. Ritz, Annals Phys. 318, 119 (2005), hep-ph/0504231. 

[10] J. Engel, M. J. Ramsey-Musolf, and U. van Kolck, Prog. Part. Nucl. Phys. 71, 21 (2013), 
1303.2371. 

[11] W. Dekens et al, JHEP 1407, 069 (2014), 1404.6082. 


22 








[12] I. B. Khriplovich, Nucl. Phys. B352, 385 (1991). 

[13] J. Engel, P. H. Frampton, and R. P. Springer, Phys. Rev. D53, 5112 (1996), nucl- 
th/9505026. 

[14] M. J. Ramsey-Musolf, Phys. Rev. Lett. 83, 3997 (1999), hep-ph/9905429, [Erratum: 
Phys. Rev. Lett.84,5681(2000)]. 

[15] A. Kurylov, G. C. McLaughlin, and M. J. Ramsey-Musolf, Phys. Rev. D63, 076007 (2001), 
hep-ph/0011185. 

[16] Y.-H. Song, R. Lazauskas, and V. Gudkov, Phys. Rev. C84, 025501 (2011), 1105.1327. 

[17] E. Mereghetti, J. de Vries, R. G. E. Timmermans, and U. van Kolck, Phys. Rev. C88, 
034001 (2013), 1305.7049. 

[18] Yu. N. Uzikov and A. A. Temerbayev, Phys. Rev. C92, 014002 (2015), 1506.08303. 

[19] P. F. Bedaque and U. van Kolck, Ann. Rev. Nucl. Part. Sci. 52, 339 (2002), nucl- 
th/0203055. 

[20] E. Epelbaum, H.-W. Hammer, and U.-G. Meifiner, Rev. Mod. Phys. 81, 1773 (2009), 
0811.1338. 

[21] E. Mereghetti and U. van Kolck, Ann. Rev. Nucl. Part. Sci. 65, 215 (2015), 1505.06272. 

[22] D. B. Kaplan and M. J. Savage, Nucl. Phys. A556, 653 (1993), [Erratum; Nucl. 
Phys.A580,679(1994)]. 

[23] S.-L. Zhu, C. M. Maekawa, B. R. Holstein, M. J. Ramsey-Musolf, and U. van Kolck, Nucl. 
Phys. A748, 435 (2005), nucl-th/0407087. 

[24] J. de Vries, E. Mereghetti, R. Timmermans, and U. van Kolck, Annals Phys. 338, 50 
(2013), 1212.0990. 

[25] E. Mereghetti, W. Hockings, and U. van Kolck, Annals Phys. 325, 2363 (2010), 1002.2391. 

[26] J. Bsaisou, U.-G. MeiBner, A. Nogga, and A. Wirzba, Annals Phys. 359, 317 (2015), 
1412.5471. 

[27] J. Dai, M. J. Savage, J. Liu, and R. P. Springer, Phys. Lett. B271, 403 (1991). 

[28] B. G. Tiburzi, Phys. Rev. D85, 054020 (2012), 1201.4852. 

[29] U.-G. Meifiner and H. Weigel, Phys. Lett. B447, 1 (1999), nucl-th/9807038. 

[30] G. ’t Hooft, Phys. Rev. Lett. 37, 8 (1976). 

[31] V. Baluni, Phys.Rev. D19, 2227 (1979). 

[32] R. Crewther, P. Di Vecchia, G. Veneziano, and E. Witten, Phys.Lett. B88, 123 (1979). 

[33] R. D. Peccei and H. R. Quinn, Phys. Rev. Lett. 38, 1440 (1977). 


23 



[34] B. Grzadkowski, M. Iskrzynski, M. Misiak, and J. Rosiek, JHEP 1010, 085 (2010), 
1008.4884. 

[35] J. de Vries et al, Phys.Rev. C84, 065501 (2011), 1109.3604. 

[36] G. Degrassi, E. Franco, S. Marchetti, and L. Silvestrini, JHEP 0511, 044 (2005), hep- 
ph/0510137. 

[37] J. Hisano, K. Tsumura, and M. J. S. Yang, Phys. Lett. B 713, 473 (2012), 1205.2212. 

[38] W. Dekens and J. de Vries, JHEP 1305, 149 (2013), 1303.3156. 

[39] J. F. Kamenik, M. Papucci, and A. Weiler, Phys. Rev. D85, 071501 (2012), 1107.3143, 
[Erratum: Phys. Rev.D88,no.3,039903(2013)]. 

[40] F. Sala, JHEP 03, 061 (2014), 1312.2589. 

[41] Y. T. Ghien, V. Girigliano, W. Dekens, J. de Vries, and E. Mereghetti, JHEP 02, 011 
(2016), 1510.00725, [JHEP02,011(2016)]. 

[42] S. Weinberg, Phys. Rev. Lett. 63, 2333 (1989). 

[43] E. Braaten, G. S. Li, and T. G. Yuan, Phys. Rev. D 42, 276 (1990). 

[44] M. J. Ramsey-Musolf and S. Su, Phys. Rept. 456, 1 (2008), hep-ph/0612057. 

[45] J. Ng and S. Tulin, Phys. Rev. D 85, 033001 (2012), 1111.0649. 

[46] E. Epelbaum, H. Krebs, and U.-G. Meifiner, Phys. Rev. Lett. 115, 122301 (2015), 
1412.4623. 

[47] V. Bernard and U.-G. Meifiner, Ann. Rev. Nucl. Part. Sci. 57, 33 (2007), hep-ph/0611231. 

[48] E. E. Jenkins and A. V. Manohar, Phys.Lett. B255, 558 (1991). 

[49] V. Bernard, N. Kaiser, J. Kambor, and U.-G. Meifiner, Nucl.Phys. B388, 315 (1992). 

[50] S. Weinberg, Physica A96, 327 (1979). 

[51] J. Gasser and H. Leutwyler, Annals Phys. 158, 142 (1984). 

[52] M. Pavon Valderrama and D. R. Phillips, Phys. Rev. Lett. 114, 082502 (2015), 1407.0437. 

[53] S. Weinberg, Phys. Lett. B251, 288 (1990). 

[54] C. Ordonez, L. Ray, and U. van Kolck, Phys. Rev. Lett. 72, 1982 (1994). 

[55] V. Bernard, N. Kaiser, and U.-G. Meifiner, Int.J.Mod.Phys. E4, 193 (1995), hep- 
ph/9501384. 

[56] E. Epelbaum, W. Glockle, and U.-G. Meifiner, Nucl.Phys. A747, 362 (2005), nucl- 
th/0405048. 

[57] N. Kaiser and U.-G. Meifiner, Nucl. Phys. A510, 759 (1990). 


24 



[58] D. R. Phillips, D. Samart, and C. Schat, Phys. Rev. Lett. 114, 062301 (2015), 1410.1157. 

[59] S.-L. Zhu, S. J. Puglia, B. R. Holstein, and M. J. Ramsey-Musolf, Phys. Rev. D63, 033006 
(2001), hep-ph/0005281. 

[60] M. Viviani et al, Phys. Rev. C89, 064004 (2014), 1403.2267. 

[61] L. Girlanda, Phys. Rev. C77, 067001 (2008), 0804.0772. 

[62] E. Epelbaum, W. Glockle, and U.-G. Meifiner, Eur. Phys. J. A19, 125 (2004), nucl- 
th/0304037. 

[63] N. Kaiser, Phys. Rev. C76, 047001 (2007), 0711.2233. 

[64] J. de Vries et al, Eur. Phys. J. A50, 108 (2014), 1404.1576. 

[65] D. R. Phillips, M. R. Schindler, and R. P. Springer, Nucl. Phys. A822, 1 (2009), 0812.2073. 

[66] H. W. Griesshammer and M. R. Schindler, Eur. Phys. J. A46, 73 (2010), 1007.0734. 

[67] B. Desplanques, J. P. Donoghue, and B. R. Holstein, Annals Phys. 124, 449 (1980). 

[68] A. Manohar and H. Georgi, Nucl. Phys. B234, 189 (1984). 

[69] J. Wasem, Phys. Rev. C85, 022501 (2012), 1108.1151. 

[70] E. G. Adelberger et al, Phys. Rev. C27, 2833 (1983). 

[71] W. G. Haxton, Phys. Rev. Lett. 46, 698 (1981). 

[72] S. A. Page et al, Phys. Rev. C35, 1119 (1987). 

[73] C. S. Wood et al. Science 275, 1759 (1997). 

[74] V. V. Flambaum and D. W. Murray, Phys. Rev. C56, 1641 (1997), nucl-th/9703050. 

[75] W. C. Haxton, C. P. Liu, and M. J. Ramsey-Musolf, Phys. Rev. Lett. 86, 5247 (2001), 
nucl-th/0101018. 

[76] E. Epelbaum, U.-G. Meifiner, W. Glockle, and C. Elster, Phys. Rev. C65, 044001 (2002), 
nucl-th/0106007. 

[77] J. de Vries, E. Mereghetti, R. G. E. Timmermans, and U. van Kolck, Phys. Rev. Lett. 
107, 091804 (2011), 1102.4068. 

[78] A. Pich and E. de Rafael, Nucl. Phys. B367, 313 (1991). 

[79] W. Hockings and U. van Kolck, Phys.Lett. B605, 273 (2005), nucl-th/0508012. 

[80] K. Ottnad, B. Kubis, U.-G. Meifiner, and F.-K. Guo, Phys.Lett. B687, 42 (2010), 
0911.3981. 

[81] J. de Vries, R. G. E. Timmermans, E. Mereghetti, and U. van Kolck, Phys. Lett. B695, 
268 (2011), 1006.2304. 


25 



[82] C.-Y. Seng, J. de Vries, E. Mereghetti, H. H. Patel, and M. Ramsey-Musolf, Phys. Lett. 
B736, 147 (2014), 1401.5366. 

[83] C. Maekawa, E. Mereghetti, J. de Vries, and U. van Kolck, Nucl.Phys. A872, 117 (2011), 
1106.6119. 

[84] J. Bsaisou et al, JEEP 1503, 104 (2015), 1411.5804. 

[85] J. Bsaisou et al, Eur.Phys.J. A49, 31 (2013), 1209.6306. 

[86] J. de Vries, E. Mereghetti, and A. Walker-Loud, Phys. Rev. C92, 045201 (2015), 
1506.06247. 

[87] A. Walker-Loud, PoS LATTICE2013, 013 (2014), 1401.8259. 

[88] S. Borsanyi et al, Science 347, 1452 (2015), 1406.4088. 

[89] M. Hoferichter, J. Ruiz de Elvira, B. Kubis, and U.-G. Meifiner, Phys. Rev. Lett. 115, 
092301 (2015), 1506.04142. 

[90] M. Pospelov, Phys.Lett. B530, 123 (2002), hep-ph/0109044. 

[91] E. Mereghetti, J. de Vries, W. Hockings, C. Maekawa, and U. van Kolck, Phys.Lett. B696, 
97 (2011), 1010.4078. 

[92] B. Borasoy, Phys.Rev. D61, 114017 (2000), hep-ph/0004011. 

[93] E. Shintani et al, Phys.Rev. D72, 014504 (2005), hep-lat/0505022. 

[94] E. Shintani, S. Aoki, and Y. Kuramashi, Phys.Rev. D78, 014503 (2008), 0803.0797. 

[95] D. O’Connell and M. J. Savage, Phys. Lett. B633, 319 (2006), hep-lat/0508009. 

[96] T. Akan, P.-K. Guo, and U.-G. Meifiner, Phys.Lett. B736, 163 (2014), 1406.2882. 

[97] P.-K. Guo and U.-G. Meifiner, JEEP 1212, 097 (2012), 1210.5887. 

[98] P. K. Guo et al, Phys. Rev. Lett. 115, 062001 (2015), 1502.02295. 

[99] A. Shindler, T. Luu, and J. de Vries, Phys. Rev. D92, 094518 (2015), 1507.02343. 

[100] M. Liischer, JEEP 1008, 071 (2010), 1006.4518. 

[101] T. Bhattacharya, V. Cirigliano, R. Gupta, E.-W. Lin, and B. Yoon, Phys. Rev. Lett. 115, 
212002 (2015), 1506.04196. 

[102] M. Pospelov and A. Ritz, Phys.Rev. D63, 073015 (2001), hep-ph/0010037. 

[103] J. Eisano, J. Y. Lee, N. Nagata, and Y. Shimizu, Phys.Rev. D85, 114044 (2012), 
1204.2653. 

[104] D. A. Demir, M. Pospelov, and A. Ritz, Phys. Rev. D67, 015007 (2003), hep-ph/0208257. 


26 



[105] S. Inoue, M. J. Ramsey-tree level, and Y. Zhang, Phys.Rev. D89, 115023 (2014), 
1403.4257. 

[106] R. B. Wiringa, V. G. J. Stoks, and R. Schiavilla, Phys. Rev. C51, 38 (1995), nucl- 
th/9408016. 

[107] V. G. J. Stoks, R. A. M. Klomp, C. P. F. Terheggen, and J. J. de Swart, Phys. Rev. C49, 
2950 (1994), nucl-th/9406039. 

[108] M. R. Schindler and R. P. Springer, Nucl. Phys. A846, 51 (2010), 0907.5358. 

[109] H. W. Griesshammer, M. R. Schindler, and R. P. Springer, Eur. Phys. J. A48, 7 (2012), 
1109.5667. 

[110] J. Vanasse and M. R. Schindler, Phys. Rev. C90, 044001 (2014), 1404.0658. 

[111] D. B. Kaplan, M. J. Savage, R. P. Springer, and M. B. Wise, Phys. Lett. B449, 1 (1999), 
nucl-th/9807081. 

[112] M. J. Savage and R. P. Springer, Nucl. Phys. A644, 235 (1998), nucl-th/9807014, [Erra¬ 
tum: Nucl. Phys.A657,457(1999)]. 

[113] D. B. Kaplan, M. J. Savage, and M. B. Wise, Phys. Lett. B424, 390 (1998), nucl- 
th/9801034. 

[114] S. Fleming, T. Mehen, and 1. W. Stewart, Nucl. Phys. A677, 313 (2000), nucl-th/9911001. 

[115] J. Carlson, R. Schiavilla, V. R. Brown, and B. F. Gibson, Phys. Rev. C65, 035502 (2002), 
nucl-th/0109084. 

[116] J. de Vries, U.-G. Meifiner, E. Epelbaum, and N. Kaiser, Eur. Phys. J. A49, 149 (2013), 
1309.4711. 

[117] D. R. Entem and R. Machleidt, Phys. Rev. C68, 041001 (2003), nucl-th/0304018. 

[118] E. Epelbaum, H. Krebs, and U.-G. MeiBner, Eur. Phys. J. A51, 53 (2015), 1412.0142. 

[119] A. Nogga, R. G. E. Timmermans, and U. van Kolck, Phys. Rev. C72, 054006 (2005), 
nucl-th/0506005. 

[120] D. R. Phillips, PoS CD12, 013 (2013), 1302.5959. 

[121] C. M. Maekawa and U. van Kolck, Phys. Lett. B478, 73 (2000), hep-ph/0006161. 

[122] M. J. Musolf and B. R. Holstein, Phys. Rev. D43, 2956 (1991). 

[123] P. D. Eversheim et al, Phys. Lett. B256, 11 (1991). 

[124] S. Kistryn et al, Phys. Rev. Lett. 58, 1616 (1987). 

[125] TRIUMF E497, A. R. Berdoz et al, Phys. Rev. Lett. 87, 272301 (2001), nucl-ex/0107014. 

[126] D. E. Driscoll and G. A. Miller, Phys. Rev. C39, 1951 (1989). 


27 



[127] D. E. Driscoll and U.-G. Meifiner, Phys. Rev. C41, 1303 (1990). 

[128] R. Schiavilla, J. Carlson, and M. W. Paris, Phys. Rev. 070, 044007 (2004), nucl- 
th/0404082. 

[129] J. F. Cavaignac, B. Vignon, and R. Wilson, Phys. Lett. B67, 148 (1977). 

[130] M. T. Gericke et al, Phys. Rev. C83, 015505 (2011). 

[131] C. Crawford, Talk given at the 3rd Workshop on the Physics of Fundamental Symmetries 
and Interactions at Low Energies and the Precision Frontier on behalf of the NPDGamma 
collaboration, PSI (2013). 

[132] R. Schiavilla, J. Carlson, and M. W. Paris, Phys. Rev. C67, 032501 (2003), nucl- 
th/0212038. 

[133] C. H. Hyun, T.-S. Park, and D.-P. Min, Phys. Lett. B516, 321 (2001), nucl-th/0101043. 

[134] C. P. Liu, Phys. Rev. C75, 065501 (2007), nucl-th/0609078. 

[135] J. de Vries et al, Phys. Lett. B747, 299 (2015), 1501.01832. 

[136] A. Cox, S. Wynchank, and C. Collie, Nuclear Physics 74, 497 (1965). 

[137] S. Rolling, E. Epelbaum, H. Krebs, and U.-G. Meifiner, Phys. Rev. C80, 045502 (2009), 
0907.3437. 

[138] Y.-H. Song, R. Lazauskas, and V. Gudkov, Phys. Rev. C86, 045502 (2012), 1207.7039. 

[139] E. Epelbaum et al, Phys. Rev. C66, 064001 (2002), nucl-th/0208023. 

[140] C. Crawford, Private communication. 

[141] J. Lang et al, Phys. Rev. Lett. 54, 170 (1985). 

[142] T. Roser and M. Simonius, Nucl. Phys. A442, 701 (1985). 

[143] S. Elhatisari et al, Nature 528, 111 (2015), 1506.03513. 

[144] W. M. Snow et al, Phys. Rev. C83, 022501 (2011). 

[145] F. J. M. Farley et al, Phys. Rev. Lett. 93, 052001 (2004), hep-ex/0307006. 

[146] JEDI-collaboration, http://collaborations.fz-juelich.de/ikp/jedi/ . 

[147] JEDI, D. Eversmann et al, Phys. Rev. Lett. 115, 094801 (2015), 1504.00635. 

[148] O. Sushkov, V. Flambaum, and 1. Khriplovich, Sov. Phys. - JETP 60:5 (1984). 

[149] O. Lebedev, K. A. Olive, M. Pospelov, and A. Ritz, Phys.Rev. D70, 016003 (2004), 
hep-ph/0402023. 

[150] C. P. Liu and R. G. E. Timmermans, Phys. Lett. B634, 488 (2006), nucl-th/0602010. 


28 



[151] Y.-H. Song, R. Lazauskas, and V. Gudkov, Phys. Rev. C83, 065503 (2011), 1104.3051. 

[152] I. Khriplovich and R. Korkin, Nucl.Phys. A665, 365 (2000), nucl-th/9904081. 

[153] C.-P. Liu and R. Timmermans, Phys.Rev. C70, 055501 (2004), nucl-th/0408060. 

[154] N. Yamanaka and E. Hiyama, Phys.Rev. C91, 054005 (2015), 1503.04446. 

[155] I. Stetcu, C. P. Liu, J. L. Friar, A. C. Hayes, and P. Navratil, Phys. Lett. B665, 168 
(2008), 0804.3815. 

[156] Y.-H. Song, R. Lazauskas, and V. Gudkov, Phys.Rev. C87, 015501 (2013), 1211.3762. 

[157] W. Dekens and D. Boer, Nucl. Phys. B889, 727 (2014), 1409.4052. 

[158] B. M. Roberts, V. A. Dzuba, and V. V. Flambaum, Ann. Rev. Nucl. Part. Sci. 65, 63 
(2015), 1412.6644. 

[159] E. Epelbaum and U.-G. MeiBner, Ann. Rev. Nucl. Part. Sci. 62, 159 (2012), 1201.2136. 

[160] J. Garlson et al, (2014), 1412.3081. 


29 



